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GOVERNMENT RIGHTS 

[0001] The U.S. government has rights in the disclosed invention pursuant to NIH Grant No. 
P41RR09784 to Stanford University. 

BACKGROUND OF THE INVENTION 

[00021 This invention relates generally to magnetic resonance imaging of an object having 
different chemical species therein, such as fat and water, and more particularly the invention 
relates to species imaging in the presence of magnetic field heterogeneity. 
[00031 The ability to image different chemical species such as fat and water using magnetic 
resonance techniques is particularly important in medical applications. For example, imaging 
and diagnosis of articular cartilage abnormalities has become increasingly important in the 
setting of an aging population where osteoarthritis is second only to cardiovascular disease as a 
cause of chronic disability. Accurate assessment of articular cartilage is also essential with the 
advent of surgical and pharmacological therapies that require advanced imaging techniques for 
initial diagnosis and management of disease progression. 

[00041 Ideal imaging of articular cartilage requires high resolution and good contrast with 
adjacent tissues; this can be markedly improved with fat suppression techniques. In addition, 
bright appearance of synovial fluid is advantageous as it provides an arthroscopic effect and "fills 
in" defects in articular cartilage, increasing the conspicuity of cartilage irregularities. Separating 
fat and water can increase the conspicuity of the both the water (for most applications) and fat 
(for special applications) with many types of pulse sequences and with both Tl and T2 weighted 
sequences. 

[00051 The difficulty in decomposing different chemical species is compounded by the 
presence of magnetic field heterogeneity. Separation of fat and water through "in-phase" and 
"out-of-phase" imaging is an approach first demonstrated by Dixon, Radiology 1984; 153: 189- 
194, and forther refined by Glover, Journal of Magnetic Resonance Imaging 1991; 1:521-530, to 
compensate for the effects of magnetic field heterogeneities. In Glover's work, a three-point 
sampling scheme that acquires spin-echo or gradient echo images with echo time (TE) 
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increments of 0, 2.2, and 4.4ms, and produce phase increments of 0, n, and 2n, when the 
frequency difference between fat and water is approximately -220Hz at 1.5T. The mathematics 
for this special case are greatly simplified and post-processing calculations are faster; however, 
these values of TE lengthen the minimum TR and would cause severe image degradation with 
SSFP imaging, for example, in the presence of typical magnetic field heterogeneities. 
Application of "Dixon" imaging to fast spin-echo (FSE) sequences has also been limited because 
the acquisition of echoes at different time shifts with respect to the spin-echo increases the 
spacing between successive refocusing pulses (echo spacing). Increasing the echo spacing 
reduces the number of echoes that can be collected in a time that maintains acceptable blurring 
from T2 decay, offsetting the scan time benefits of FSE. A fat-water separation method that 
permitted shorter time increments would reduce the time between refocusing pulses and be 
beneficial to fast spin-echo imaging. 

[00061 SSFP is a rapid gradient echo imaging technique with renewed interest in recent 
years, owing to widespread availability of high speed gradient systems. SSFP has superior signal 
to noise ratio (SNR) compared to other gradient echo techniques and has excellent contrast 
behavior that has mixed dependence on Tl and T2. In particular, synovial fluid appears bright on 
SSFP images owing to its long T2. The major limitation of SSFP is severe image degradation 
caused by local magnetic field heterogeneities if the repetition time (TR) is long. 



SUMMARY OF THE INVENTION 

[0007] The present invention utilizes a multi-point chemical species separation process which 
is compatible with a rapid gradient echo imaging technique, such as SSFP, fast spin echo, echo 
planar imaging, spin echo, spiral imaging, and other similar pulse sequences. An iterative least 
squares fitting algorithm is utilized to combine signals at different echo times using an assumed 
initial value of field heterogeneity. A calculated value of field heterogeneity error is then 
obtained from the combined signals, and a new value of combined signals is obtained using the 
initially assumed value of field heterogeneity and the calculated error. An updated value of error 
in the field heterogeneity estimate is calculated and the process is repeated until an acceptable 
error value is realized. 

[0008] The linear combination technique facilitates the estimation of water and fat images 
from a minimum of three images acquired at different, but arbitrarily spaced, echo times. If the 
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field heterogeneity map, i^, is known, a minimum of two images acquired at different echo times 
are sufficient to estimate water and fat images. 

[00091 The method of imaging different chemical species, such as water and fat, can be 
implemented with a single coil or with multiple coils. A field map of Bo heterogeneity is first 
calculated for each coil, and then the field maps are combined by weighting the contributions 
fi-om each coil, such as by the square of the magnitude of the image contributed fi-om each coil. 
Using the combined field map, images firom each of the coils are recalculated, and the 
recalculated images are then combined, such as by taking the square root of the sum of the 
squares of the individual coil image signals. The combined field map can be smoothed by low 
pass filtering before the images are recalculated in order to improve noise performance. 
[00101 The invention and objects and features thereof will be more readily apparent firom the 
following detailed description and appended claims when taken with the drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[001 1] Fig. 1 A is a flow diagram illustrating a magnetic resonance imaging using a single 
coil in accordance with one embodiment of the invention. 

[00121 Fig. 1 B is a flow diagram of another method of magnetic resonance imaging using a 
plurality of coils in accordance with another embodiment of the invention. 
[00131 Figs. 2A-2D are images of knees constixicted using the invention and using spoiled 
gradient echo (SPGR) imaging for comparison at 1.5T. 

[00141 Figs. 3A-3C are images of knees constructed using the invention at 3.0T. 
[00151 Figs. 4A-4C are images of ankles constructed using tiie invention at 1 .5T. 
[001 6] Figs. 5 A-5D are images of knees constiiicted using tiie invention and using spoiled 
echo (SPGR) for comparison. 

[00171 Fig. 6 is cardiac images using the invention and acquired throughout the cardiac cycle 
including source images and calculated water and fat images using the invention. 
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DETAILED DESCRIPTION OF ILLUSTRATED EMBODIMENTS 

[0018] The present invention provides a new method for implementing fat-water separation 
utilizing a novel iterative least-squares method that reconstructs data acquired at short echo time 
increments, yielding images with high SNR and uniform separation of fat signal from water 
signal. The algorithm summarized in Fig. 1 extends naturally to multi-coil reconstruction with 
minimal additional complexity as summarized in Fig. 2. Single and multi-coil decompositions 
derived from images obtained at both 1 .5T and 3.0T are shown. Examples in the knee, ankle, 
pelvis and heart are shown. An analysis of the noise performance of this method is provided and 
methods for improving noise performance through field map smoothing are discussed. 
[00191 As noted above, the use of short echo times is necessary with SSFP to maintain short 
TR, thereby preventing image degradation from field heterogeneities. Three point methods 
described previously by Glover describe the special case of phase shifts of 0, n, and 2;r which 
correspond to echo tune increments of 0, 2.2ms and 4.4ms at 1.5T and 0. 1.1ms and 2.2ms at 
3 0T. Unfortunately, such echo time increments cause significant lengthening TR. When using 
SSFP, increases in TR can lead to severe image degradation from banding artifacts caused by off- 
resonance field heterogeneities and chemical shift. Decomposition of water and fat using 
"Dixon" approaches requires shorter increments in echo time. 

[00201 In accordance with the invention, an iterative linear least-squares approach is 
formulated and a generalized algorithm with arbitrary echo times and multiple chemical species 
is described below. The invention is extended to multi-coil applications, and an algorithm 
apphes the invention to decomposition of each chemical species. 

[00211 Consider the signal in an image from a pixel containing M species each with chemical 
shift 4/; (Hz) located at position r, at an echo time t. 



I M 



(1) 



where pj is the signal from the/ species and is. in general, a complex term with its own 



magnitude, \pj 1 and phase, and m is the local magnetic field heterogeneity (Hz). If 
measurements are made at discrete echo times tn («=1, ••. ^0. then. 
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" - ^J2^ir),„ (2) 



representing the signal in a pixel located at position r at echo time f„. Equation (2) contains M 
complex unknowns (pj ...M)) and one scalar unknown {y/), for a total of 2M+1 unknowns. 
Each image contains a real and imaginary component, constituting two measurements per time 
point t„. In general, M+1 or more images are required to determine the system and separate each 
chemical species. For example, with fat and water only (M=2), at least three or more images are 
required to decompose fat from water. If the relative phase between all pj were known, then only 
M images (=2 for fat and water, only) would be required for decomposition. 
(00221 If an initial estimate of the field map Wo(r) is known, then equation (2) can be 
rewritten. 



^ ^ J^'^j'- (3) 



and is a linear system of complex equations that can be split into real (f f ) and imaginary 



parts, 



5- = + is'„ = f:(p;c,„ - p'jdj„ )+ + p>,„ ) (4) 

7=1 



where p] and p'j are the real and imaginary components of the/ species, Cj„ = cos(2;rAr/„ ) 
and dj„ = sm(27ti^jt„ ). Equation (4) forms a set of linear equations that is amenable to linear 
least squares fitting to decompose each chemical species. For n=\, .... N, equation (4) can be 
written in matrix format, 

(5) 

S = Ap 

where S = [5f ... ... S',Y,p = W Pi P'^- P'-^ '"^^ 

matrix A is provided below for M species: 



STFUP152/S02-130 



Page 5 



A = 









"21 












^22 


^22 




^A/2 


"a/2 














J 


dn 






^21 








dn 




^22 


^22 




^Af2 



















(A.1) 



*2W 



-2Ar 



'MAT 



where the rows 1 to are used to calculate the real components of the signal, and rows N+1 to 
2N are used to calculate the imaginary components of the signal. 

[0023] Using well-described least squares fitting approaches for linear systems of equations, 
it can be shown that the least-squares estimate of p is. 



[00241 From equation (6), initial estimates of the water and fat images can be determined 
fromp = [y5f p( M M - p'mJ ■ 

[0025] The initial estimate of the field map {y/Jis further refined by defining error terms: 
= ^„ + Av^ , = y3; + Apf , and /?; = p'j + Ap'j . Liserting these expressions into equation 
(2) it is shown in the appendix that for small Ay/ , Apj and Apj , that, 

where y = [A<^ Apf Ap[ Ap* Api ... Apl Ap^, f , and 1 and B are both defined in 

the appendix. For «=1 equation (7) is a linear system of equations, and similar to above, 

estimates of y can be calculated as, 

y = (B^Br'B^i (8)- 

[00261 In the special case of a system with only water and fat ( A/^ chemical shift) and the 

receive/transmit frequency of the scanner is set to water, matrices A and B become. 
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A = 



1 
1 

1 
0 
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0 
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1 
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0 



d^ 



and 



B = 



^1 
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0 
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d^ I 
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d^ 



d^ 



'H . 



where = cos(2;rA/^rJ, ^> = sin(2M/^r J, -^ntX' pl~ P]d^ - p\c^ and 
= 2nt^ (pi -^PfC,- p'fd„ ) are the matrix elements. 

[0027] Using the above equations, the following algorithm summarizes the method used to 
determine the least-squares estimates of different chemical species, such as water images, fat 
images, and silicone images for each pixel using a single coil acquisition: 

a) estimate water and fat assuming with equation (6), and initial guess for field map, y/o^ A 
useful initial guess is zero (Hz); 

b) calculate error to field map, ^t^^ using equation (8); 

c) recalculate i// = if/^+Ai// ; 

d) recalculate S = [i, ... s^Y with the new estimate ofy/ using equation (6); 

e) repeat steps b) through d) until zl^^is small (e.g., <1 Hz). 
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[0028] A multi-coil acquisition with P elements collects P independent images, all with a 
different relative phase offset. By using the above algorithm, P separate water and fat images can 
be generated, as well as P field heterogeneity maps. The field maps can then be combined by 
weighting the contribution from each coil by the square of the magnitude of the image 
contributed by that coil. Specifically, for each pixel at position r, the combined field map can be 
calculated as, 



(11). 



[0029] Using the combined field map, S = Se^'^'"^'^^'^^ is then recalculated and estimates of the 
water and fat images are calculated for each coil from equation (6). Finally, the M water and fat 
images are combined using the square root of the sum of the square, a commonly used multi-coil 
reconstruction. 

[0030] A summary of the final fat-water decomposition is provided below, taking into 
account the use of multi-coil acquisitions, as well as smoothing of the final field map: 

For each coil: 

for each pixel; 

a) estimate water and fat assuming initial guess for field map, y/o , from equation (6); 

b) calculate error to field map, Ay/^ from equation (8); 

c) recalculate t// = Ay/ ; 

d) recalculate S = [i, ^2 ^nV ^^^h the new estimate ofy/ ; and 

e) repeat steps b) through d) until Ay/ is small. 

For multi-coil acquisitions, combine field map with equation (11). 
Smooth field map with low-pass filter. 
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Recalculate final estimate of each water and fat image with equation (6) using the filtered 
combined field map. 

Multi-coil combination of water (fat) images using the square root of the smn of the squares for 
all signals. 

[0031] A three dimensional (3-D) SSFP imaging sequence was used on a 1 .5T and a 3.0T GE 
Signa scanner to acquire source, calculated water and calculated fat images from knees (Figs. 
2A-2D, Figs. 3A-3C), ankles (Figs. 4A-4C, Figs. 5 A-5D) and abdomens (Fig. 6) from normal 
volunteers. 

[0032] A product extremity coil was used for the knees and ankles and the body coil was 
used for abdominal imaging. Pulse sequence parameters included: BW=125kHz, Nx=256, 
Ny=192, NSA=1. For the knee and ankle: F0V-16cm, slice=1.5mm, Nz=64, TR=6.2, TE=(2.02, 
2.69, 3.35, 4.02ms), total scan time of 5:02 min. Pelvis: FOV=32cm, slice=5mm, N2=22, 
TR=5.4ms, TE=(1.58, 2.24, 2.91,3.58ms), total scan time of 1:30 min. 

[0033] For comparison, 3D spoiled gradient echo (SPGR) images with fat saturation were 
obtained in the knees and ankles (Figs. 2D, 5D). Matrix size, field of view and slice thickness 
were all the same. Other parameters included: TR=50ms, TE=5ms (full echo), tip angle = 40*^, 
and BW = ± 16kHz. These parameters are based on established reports using fat saturated SPGR 
imaging of articular cartilage. 

[0034] Figs. 6A-6I are retrospectively ECG gated cardiac CINE SSFP images acquired at 
end-diastole, mid-systolic, and end-systoUc, respectively, for source images (Figs. 6A-6C), 
calculated water images (Figs. 6D-6F), and calculated fat images (Figs. 60-61). 

[0035] An off-line reconstruction program written in Matlab 6.0 (Mathworks, Mountain 
View, CA) was used to perform fast Fourier transform reconstruction of all images. Following 
reconstruction of complex (magnitude and phase) images, estimation of water images and fat 
images based on the iterative least-squares algorithm was performed. 

[0036] Multi-point techniques have been combined with short TR and TE sequences such as 
SSFP to obtain excellent fat-water separation, using iterative least-squares fitting approaches that 
allow the use of short echo time increments. This invention is particularly attractive because it 
does not require evenly-spaced echo time increments. In addition, phase-unwrapping algorithms 
are not necessary. While the invention works well with short TE/short TR sequences, the method 
works with other sequences in "Dixon" imaging. 
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[0037] Incorporated herein by reference for all purposes and attached hereto as Appendix B 
is a copy of a paper submitted to Magnetic Resonance in Medicine which has been accepted for 
publication. The paper is entitled "Multi-coil "Dixon" Chemical Species Separation with an 
Iterative Least Squares Estimation Method", by Scott B. Reeder, Zhifei Wen, Huanzhou Yu, 
Angel R. Pineda, Garry E. Gold, Michael Marki, and Norbert J. Pelc. 

[0038] While the invention has been described with reference to specific embodiments, the 
description is illustrative of the invention and is not to be construed as limiting the invention. For 
example, rather than fiUering the combined field map, the source images can be smoothed by 
fihering before calculating the field map. Then with this field map, go back to the original 
unfilitered source images and calculate the water and fat images. Thus, various modifcations and 
applications may occur to those skilled in the art without departing from the true scope and spirit 
of the invention as defined by the appended claims. 
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